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According to Lipatov, high orders of perturbation theory are determined 
by the saddle-point configurations (instantons) of the corresponding func- 
tional integrals. According to 't Hooft, some individual large diagrams, 
renormalons, are also significant and they are not contained in the Lipatov 
contribution. The history of the conception of renormalons is presented, and 
the arguments in favor of and against their significance are discussed. The 
analytic properties of the Borel transforms of functional integrals, Green 
functions, vertex parts, and scaling functions are investigated in the case of 
ip 4 theory. Their analyticity in a complex plane with a cut from the first 
instanton singularity to infinity (the Le Guillou- Zinn- Justin hypothesis) is 
proved. It rules out the existence of the renormalon singularities pointed 
out by 't Hooft and demonstrates the nonconstructiveness of the conception 
of renormalons as a whole. The results can be interpreted as an indication 
of the internal consistency of ip 4 theory. 
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1. INTRODUCTION 



Many problems in theoretical physics can be reduced to a calculation of functional 
integrals of the type 

I = J D<pexp(-S {<p} - gS int {(p}), (1) 

whose expansion in the coupling constant g gives an ordinary perturbation theory. In 1977 
Lipatov pQ proposed a method for calculating the high-order expansion coefficients of the 
integrals (1) on the basis of the following simple idea. If the function F(g) is expanded into 
a series 

oo 

F(g) = £ F N g N , 

N=0 

the iVth expansion coefficient can be calculated from the formula 

f dg F{g) 

/ 2™^' (2) 

c 

where the contour C encloses the point g = in the complex plane. Taking the integral 
(1) as F(g), we obtain 



l N-l 



27r? / dg I D V ex P(- s o{f} ~ gSintif} - Nlng), (3) 

and the appearance of an exponential with a large exponent indicates that the saddle-point 
method can be applicable at large N. Lipatov's idea is to seek the saddle point in (3) with 
respect to g and tp simultaneously: such saddle-point exists for all the cases of interest and 
is realized on a spatially localized function, which has been termed an instanton. Moreover, 
the conditions for applicability of the saddle-point method are satisfied at large N. 

The Lipatov technique, which was originally applied to scalar theories, such as if 4 theory 

SoW} + gStofr} = J d d x |i(V^) 2 + imV + \gV 4 } , (4) 

was subsequently generalized to vector fields [2] , scalar electrodynamics [HI H] , Yang-Mills 
fields [5], fermion fields [6], etc. (see the collection of articles in Ref. p]). The ultimate goal 
was to apply it to theories of practical interest, viz., quantum electrodynamics 0, E] and 
quantum chromodynamics (QCD) [TQ1 [H] - As was pointed out already in Lipatov's first 
paper [T], the knowledge of the first few coefficients and their asymptotic behavior permits 
approximate reconstruction of the Gell-Mann-Low function, opening up a direct route to 
the solutions of the problem of confinement and electrodynamics at short distances. 

However, a conception which raised some doubts regarding the Lipatov technique was 
conceived already in 1977. It was initiated by a paper by Lautrup [T2] . which contained 
the following curious observation. The typical result of calculations based on the Lipatov 
technique has a functional form 



in 



ca N T(N + b)^ ca N N b - l N\ (5) 
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Figure 1: Example of a diagram for quantum electrodynamics, which makes a contribution 
~iV! to the iVth order of perturbation theory |12j . 

and the natural interpretation of it is that there is a factorially large number of diagrams 
of the same order (ag) N . However, in the general case such an interpretation is incorrect, 
since there are examples of individual iVth-order diagrams having a value ~iV!. The latter 
are diagrams (Fig. 1) which contain long chains of "bubbles" . Such factorial contributions 
of individual diagrams were termed renormalons, since they appear only in renormalizable 
theories. Lautrup's example (Fig. 1) was related to quantum electrodynamics, but similar 
diagrams exist in QCD and four-dimensional ip 4 theory. 

Strictly speaking, nothing followed from Lautrup's observation: the Lipatov technique is 
based on a formal calculation of the functional integral (3) and does not rely in any way on 
a statistical analysis of diagrams. It is natural to expect that the renormalon contributions 

1 In a broader sense, a renormalizable theory is one in which the divergences are eliminated by renor- 
malizing a finite number of parameters. According to more precise terminology, such theories are subdi- 
vided into super-renormalizable (renormalizable "with a surplus" ) and renormalizable in the narrow sense 
(marginally renormalizable); the latter, which gave their name to renormalons, are characterized by the 
logarithmic situation, which is needed for the appearance of factorial contributions (see below). 
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have already been taken into account in the instanton result (5). In fact, no far-reaching 
claims were made in Ref. [12J or in the relevant publications appearing shortly thereafter 

However, the tone of the publications subsequently changed dramatically. The reason 
was 't Hooft's lecture [15J, which was delivered in the same year, 1977. The term "renor- 
malon" was used in it for the first time, and it was asserted that renormalons are not 
contained in the instanton contribution (5). The authors of the subsequent publications 
[16J [30] considered 't Hooft's opinion to be self-evident and did not trouble themselves 
with argumentation. 

A convenient language for discussion, viz., the analytic properties of Borel transforms, 
was introduced in 't Hooft's lecture. The Borel transformation 



OO p 00 



F{9) = E F N9 N = E TFf / dxx»e-*g = / dxe~* £ ^(gxf 

N=0 N=0 n n ^=0 ' 



which factorially improves the convergence, is widely used in the theory of divergent series 
[30J. It is convenient to rewrite it in the form 

oo 

F(g)= Jdxe- x B(gx), (6) 
o 

OO 771 

N=0 iv • 

by introducing the Borel transform B(z) of the function F(g). The Borel transform for a 
function with the expansion coefficients (5) 

Bj(z) = E ca N N b - 1 z N ~ (1 - az)-\ za -> 1 (8) 

N 

has a singularity at the point z = 1/a. 

't Hooft arrived at this result in a different way, without reference to the Lipatov 
technique. Rewriting the integral (1) and the definition of the Borel transform (6) in the 
form 

/= D<pexp(-S{<p}/g), (9) 



>3C 

F{g) = J dxe- x/9 B{x), (10) 
o 

which can be accomplished by means of the replacements ip — > if/ \fg and x — > x/g in (4) 
and (6),0 yields the Borel transform of the integral (9): 

B ,W = /iV(»-sM)= / i^f, (ID 

z=SW} 



2 't Hooft omitted the factors of the form g n , since integrals of the type (1) usually appear in the form 
of a ratio and such factors cancel out. 
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Figure 2: More general class of renormalon diagrams 



where the latter integration is carried out over the hypersurface z = S{ip}. If an instanton 
(f c (x), i-e., a classical solution with a finite action, exists for the integral (9), then SS{<p c } = 
and the partial derivatives dS/dtpi with respect to all the variables y?j comprising Dip 
vanish; therefore, VS{<p c } = and the Borel transform (11) has a singularity at the point 



which coincides with 1/a. In addition, there are singularities at the points mS{ip c }, which 
correspond to solutions containing m infinitely distant instantons. If it is assumed that 
the singularity (12) is closest to the origin of coordinates, the result (5) of the Lipatov 
technique is reproduced. However, 't Hooft allowed the existence of singularities differing 
from those of the instanton type: in this case the asymptotics of the expansion coefficients 
can be specified by the singularity which is closest to the origin of coordinates. 

't Hooft regarded renormalons as a possible mechanism for the appearance of the new 
singularities. Let us take an arbitrary diagram for quantum electrodynamics and single out 
the line of a virtual photon with the momentum k (or an interaction line in ip A theory) in 
it (Fig. 2, a): it corresponds to integration over a region of large momenta of the type 



where m is an integer. If we assume that all the renormalizations have been performed, 
the integral converges and m > 3. Inserting a chain of iV "bubbles" into the photon line 



z = S{ip c }, 



(12) 




k>ko 
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(Fig. 2,b), we obtain the integral El 



oo , v N 

d 4 kk- 2m ((3 lnk 2 ) N ~(3» = f dtt N e^ m - 2 ^ ~ I ) M (13) 

fe>fc x 7 

Borel summation of a sequence of such diagrams gives singularities at the points 

" m = 3,4,5,.... (14) 



The constant /3q is the first nonvanishing expansion coefficient of the Gell-Mann-Low func- 
tion (Sec. 2), and with consideration of the sign relationships (S{(p c } < 0, (3o > 0) 't Hooft 
arrived at the picture of singularities for </> 4 theory shown in Fig. 3, a. 

It is not difficult to see that 't Hooft 's arguments regarding renormalons leave some 
fundamental questions unanswered: 

Why can significance be attached to individual sequences of diagrams, which make up 
an infinitesimal fraction of their total number when N — ► oo? 

How do we know that the renormalons have not already been taken into account in the 
instanton contribution (5)? 

However, the general setting the question on the possible contributions of a noninstanton 
nature to the asymptotics of the expansion coefficients has a sense: it uncovers a gap in 
the mathematical foundation of the Lipatov approach. Indeed, let the function f(x) have 
a sharp maximum at the point xq and a slow tail at large values of x (Fig. 4,), so that 
the contributions to the integral / f(x) dx from the vicinity of the maximum and from the 
tail region are comparable. An investigation of the integral for a saddle point discloses a 
maximum at xq and (provided it is sufficiently sharp) the formal applicability of the saddle- 
point method; however, a calculation of the integral in the saddle-point approximation will 
be erroneous, since the contribution of the tail will be lost. If such tails are present in the 
integral (3), the Lipatov technique can be incorrect. @ 

The essential lack of nonsaddle-point methods for calculating functional integrals makes 
it impossible to straightforwardly investigate the contribution of possible tails. But are 
there constructive arguments pointing to their existence? In principle, such arguments ex- 
ist, but they have a fairly intuitive and ambiguous character and do not hold up to criticism 
when they are closely examined (Sec. 2). As a result, the conception of renormalons has 
been in a dialectic equilibrium, i.e., it has not been proved or refuted. This uncertainty 
has caused the interest in high orders of perturbation theory to drop sharply and Lipatov's 



3 In quantum electrodynamics and QCD a polarization loop gives the factor fc 2 lnfc 2 , and a photon 
(gluon) propagator gives 1/fc 2 ; in four-dimensional (p 4 theory a closed loop corresponds to lnfc 2 , and an 
interaction line corresponds to a constant. In all cases a chain of N "bubbles" corresponds to (lnA: 2 )^. 

4 The validity of the saddle-point method can be substantiated for convergent finite- multiplicity integrals 
of functions exp[\F(x)] in the limit A — > oo (Ref. [SI]). The integral (3) can be brought into the form 
indicated, but, generally speaking, it contains both ultraviolet divergences and divergences associated with 
an infinite number of integrations. The ratio of two integrals of type (1) must be finite (after the appropriate 
renormalizations) , but each of them taken individually can be divergent. 



6 



o — o 



S{cp c } 

— o— 



mstantons 



-x — x — x — x- 

renormalons 



s i<P c } 



© 



Figure 3: a) Picture of singularities for ip 4 theory given by 't Hooft pj) 
region according to the results of the present work. 



b) Analyticity 
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Figure 4: Example of a function for which the saddle-point method is formally applicable, 
but gives a wrong result. 

program [lj to remain uncompleted. For example, a preliminary result for quantum elec- 
trodynamics was obtained back in 1978 [§], but the parameters b and c in the asymptotics 
In — ca N T(N/2 + b) have not yet been calculated. Moreover, the first result for QCD 
appeared in 1991 (Ref. [TO] ) and was recently revised [ill], but it is still unsatisfactory 
(Sec. 4), although the foundation for such calculations was completely ready in 1980 [5, 6J. 
Finally, the attempts to reconstruct the Gell-Mann-Low function have been restricted to 
ip A theory [331 EH EE] . 

A reawakening of interest in asymptotic estimates has recently been observed, but it 
has been confined almost exclusively to the renormalon doctrine [ZJJ [30]- In particular, 
it is generally accepted (see Zakharov's review [21]) that renormalons determine the per- 
turbation asymptotics in QCD. However, the work within the renormalon approach has 
already raised some doubts: the summation of larger sequences of diagrams leads to dra- 
matic renormalization of the renormalon contribution and renders the common coefficient 
in front of them totally indefinite [30J; in fact, it is impossible to state that it does not van- 
ish. On the other hand, the use of the Lipatov technique has provided significant progress 
in the theory of disordered systems [36] and in the theory of turbulence |37j . 

This paper presents a detailed discussion of the existing arguments in favor of renor- 
malons, which are shown to be unsound (Sec. 2). The analytic properties of Borel trans- 
forms are investigated in the example of y? 4 theory (Sec. 3), and their analyticity in a 
complex plane with a cut from the first instanton singularity to infinity is demonstrated 
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(Fig. 3,b). It rules out the existence of the renormalon singularities indicated by 't Hooft 
(Fig. 3, a) and demonstrates the nonconstructiveness of the conception of renormalons as a 
whole. 

An hypothesis that Borel transforms have the analytic properties indicated was ad- 
vanced by Le Guillou and Zinn- Justin [3H] and underlies one of the most efficient methods 
for summing perturbation series, which is based on the use of conformal transformation. 
The results obtained below provide the mathematical foundation of this method. 



Let us discuss the arguments in the literature that point to the existence of noninstanton 
contributions in the integral (3). 

1. There have been numerous semi-intuitive assertions which reduce to the notion that 
instantons do not exhaust all of physics. This thesis is correct as long as it is understood 
correctly, but in the present case it is not relevant. 

Historically, instantons first appeared when the saddle-point approximation was em- 
ployed in the original integral (1). It was substantiated only in a narrow region of parame- 
ters, and thus instantons did not, in fact, exhaust all of physics. In the Lipatov technique 
the situation changed, because the saddle-point approximation is used not in the integral 
(1), but in the expression (3) for the expansion coefficients. Since only large values of 
N are considered, only a limited role is assigned to instantons from the onset; however, 
the saddle-point method is now always applicable, and there is a basis to assume that 
everything is determined by instantons. 

Let us illustrate the foregoing statements in the example of the Schrodinger equation 
with a random potential V(x): 



At large negative values of E its eigenfunctions (Fig. 5) are localized on the infrequent 
fluctuations of the random potential (a and b), at large positive values of E they are similar 
to plane waves (d and e), and in the vicinity of the bare spectrum edge at E = they 
have a highly broken, fractal character (c). The problem of investigating Eq. (15) can be 
reformulated in the language of an effective field theory, viz., (f 4 theory with the "incorrect" 
sign for g (Refs. and |39J). In this case the typical wave functions of localized states 
are described by instantons. The changes in the situation observed as E increases can be 
described in the following manner in terms of instantons: at first instantons have a small 
radius and a sparse distribution, i.e., they form an ideal gas (a); then the radius of the 
instantons increases, and their density rises, i.e., interactions between them appear (b); 

5 Of course, instantons exist only in a part of the region of parameters, but this is not a restriction in 
the Lipatov technique: the values of a, b, and c in the asymptotics (5) are calculated exactly, and they 
allow analytic continuation as functions of the physical parameters. 



2. PROS AND CONS 




(15) 
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Figure 5: On the left — density of states v(E) in the presence (solid curve) and absence 
(dashed curve) of a random potential. On the right — schematic representations of the 
eigenfunctions of Eq. (15). 
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then condensation of the instantons takes place (c), and an instanton crystal forms (d and 
e). Only the case in Fig. 5, a corresponds to applicability of the saddle-point method in an 
integral of the type (1), and thus the standard instanton approximation is very poor. 

Let us examine this situation from the standpoint of perturbation theory with respect 
to the random potential V(x). An ideal instanton crystal (e) corresponds to a plane wave, 

1. e., the zeroth order of perturbation theory. In a nonideal crystal (d) the higher orders 
have an increased role; in the vicinity of the bare spectrum edge (c) all the diagrams are of 
the same order of magnitude, so that the high and low orders of perturbation theory are 
equally significant. In the region of localized states (a and b) the dominant role shifts to 
the high orders: these states are not manifested in any finite order of perturbation theory, 
and discarding the low-order contributions does not influence their properties in any way. 
We see that Lipatov's conception (high orders are determined by instantons) fits excellently 
into the existing physical picture. 

Thus, the status of instantons in the integral (1) and the integral (3) differs significantly. 
In our opinion, this accounts for the position taken by 't Hooft, since he is a classic on 
instantons [ID] specifically in the original integral (1). 

2. Relationship to the logarithmic situation [15J. Renormalons exist only in 
renormalizable theories, but not in super-renormalizable theories. If a theory is super- 
renormalizable, an upper bound of the type a N g N can be obtained for the contribution of 
an individual diagram, and the appearance of the factor N\ in the asymptotics (5) can be 
associated only with a factorially large number of diagrams. Renormalons and, thus, a new 
mechanism for the appearance of factorial contributions appear in renormalizable theories. 
It can be expected that this mechanism is associated with the formation of the tails in the 
integral (3) and is not taken into account in the Lipatov technique. 

In this argument everything except the last conclusion is correct. We can illustrate this 
in the case of ip A theory, which is renormalizable for d = 4 and super-renormalizable for 
d < 4. Among the large set of integrations concealed in the symbol Dip in the integral (3), 
we can single out one for which the limit d — > 4 is associated with qualitative changes: it 
is the integration over the instanton radius R (Fig. 6). For d significantly smaller than 4 
(for example, d = 3), the integrand exp(— S{(f}) has a sharp maximum as a function of R 
and allows saddle-point integration; when d = 4 — e the maximum becomes gently sloping, 
and when d = 4 the instanton action S{ip c } does not depend on R at all. In the latter 
case the integral diverges, leading to the logarithmic situation. We see (see the curve for 
d = 4 — e in Fig. 6) that the "activation" of renormalon contributions is, in fact, related to 
the appearance of slow tails in the integral (3), but these tails are taken into account in 
the Lipatov technique [36J. 

One can be puzzled by a following question. If the Lipatov technique is based on a 
saddle-point method, then how can it cover the definitely nonsaddle-point situation for 
d = 4 — e? The fact is that the saddle point in a functional integral practically never 
reduces to a simple maximum achieved at a single point: the maximum is degenerate 
in a certain space of finite dimensionality. Accordingly, a finite number of integrations 
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should be performed exactly, rather than in the saddle-point approximation. However, 
if the integration is performed exactly over a certain variable (for example, R), it is of 
no significance whether the degeneracy is exact (d = 4) or approximate (d = 4 — e). 
Nevertheless, in the latter case technical difficulties arise, and the corresponding methods 
(constrained instantons [S1H2]) have been poorly developed hitherto |36j . 

It is thus clear that even in cases where slow tails actually appear in the integral (3), the 
Lipatov procedure is sufficiently flexible and contains broad possibilities for dealing with 
them. 

3. The limit n — > oo. There is an opinion that the significance of renormalon contri- 
butions can easily be proved by treating the n-component (p A theory in the limit n — > oo 
(and the analogous models in QCD and electrodynamics) [16]: the factor n corresponds to 
a closed loop, and renormalon graphs containing the maximum possible number of loops 
are singled out by the large parameter n. Although diagrams of the same order, but with 
a smaller number of loops, can make comparable contributions at large N due to the com- 
binatorial factors, they have a slower dependence on n; therefore, the renormalons cannot 
be cancelled identically. 

This argument is valid in any finite order of 1/n. However, a detailed investigation of 
the structure of the 1 /n-expansion [TBI US] reveals the presence of numerous cancellations, 
and although the situation cannot be totally elucidated, the question is not resolved on the 
level of simple arguments of the type indicated. 

It is not difficult to identify the crux of the problem here. As an example, let us consider 
the self-energy E(p, m) of <£> 4 theory; it is clear from a diagrammatic analysis for m = and 
values of the momentum p close to the ultraviolet cutoff A that the (N + l)th expansion 
coefficient for E(p, 0) — E(0, 0) has the form of a polynomial in n 

p 2 {A N (N)n N + A N _ 1 (N)n N ~ 1 + ... + A 1 (N)n + A (N)}, (16) 

in which the coefficient A^(N) is specified by renormalon graphs: 

A W (A) = const- (-^)V (17) 

If it is assumed that the renormalon graphs are contained in the instanton contribution, the 
expression (16) should transform into the Lipatov asymptotics at large N [see Eq. (130) 
in Ref. 03] for M = 1 and p « A]: 

/ 1 \ N 
p^a(3 n N^/ 2 (-—-) N\ 
\ 8tc 2 J 

o 

where a,/3 ~ 1, and K\(x) is the McDonald function. It is easily seen that an equal- 
ity between (16) and (18) is impossible when n — > oo. This is a manifestation of the 
"noncancelability" of renormalons. 



r 



n + 2 
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However, the usual condition for applicability of the Lipatov technique, N ^> 1 at large 
n is, generally speaking, replaced by a more rigid condition, for example, N ^> n, and n 
then has a bound of the type 

n<n (N), (19) 

which precludes going to the limit n — > oo. If it is taken into account that the Lipatov 
asymptotics has limited accuracy (~l/iV in relative units), the correct formulation of the 
question is as follows. Can we construct an interpolation polynomial of type (16) with the 
high-order coefficient (17) which would approximate the function (18) within an assigned 
accuracy in the interval < n < n (N), where no(N) — > oo as N — > oo? The answer to 
this question is positive (see the Appendix); therefore, the assumption that the renormalon 
graphs are contained in the instanton contribution does not lead to contradictions. 



4. Relationship to a Landau pole [16], I19j . It is easy to see (Fig. 2,c) that the 
summation of a sequence of renormalon diagrams corresponds to "dressing" the interaction. 
The relationship between the renormalized charge g and the bare charge go is then given 
by the familiar expression [4*41 H5j H6] 



90 ~ l-/? sln(A7m 2 )' (20) 



which contains a pole at the point 



A 2 = m 2 e 1/l3 ° 9 . (21) 



Under the literal understanding of this pole in the spirit of the early papers by Landau and 
Pomeranchuk |17j . a simple physical interpretation can be given to renormalon singularities 
(Refs. PS] and PS]).@ 

The dependence of the perturbation series on the cutoff parameter A has the structure 

c_ x A 2 + c In A 2 + ciA" 2 + c 2 A" 4 + . . . + c n A" 2n + . . . . (22) 

The first two terms are eliminated by a renormalization procedure, and the rest terms, 
in principle, remain, but vanish in the limit A — > oo. Because of the pole in (20), values of 
A greater than A c are inaccessible in principle, and unremovable uncertainties of the type 

A; 2n = mV n//3 ° 9 (23) 

appear in the theory. Similar uncertainties are generated by renormalon singularities, 
whose existence on the positive semiaxis leads to ambiguity in the choice of the integration 
contour in the Borel integral (6). The contour can be drawn to the right or the left of the 
nth singularity (Fig. 7), producing an uncertainty in reconstructing the function from its 



6 It is assumed below that (3q > 0. For asymptotically free theories, in which (3q < 0, similar arguments 
are valid in regard to so-called infrared renormalons. The latter are obtained from integrals of the type 
(13) with ttl — 1,0, — 1, —2, ... in the region of small momenta. 
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Figure 7: Ambiguity of the choice of the integration path in the Borel integral (6) in the 
presence of renormalon singularities. 

Borel transform: 



which, with consideration of the equality z n = n/(3 coincides with (23). 

Of course, the literal interpretation of the Landau pole seems archaic, but after some 
modification of the argument presented, a real meaning can be assigned to it. It is well 
known |46j, that the dependence of the charge g on the distance scale A -1 is given by the 
equation 



whose solution depends drastically on the behavior of the Gell- Mann-Low function f3(g). 
The pole in (20) is eliminated, if (3(g) changes sign or behaves as g a with a < 1 at large 
g. If, on the other hand, (3(g) is positive and increases as g a with a > 1 when g — > oo, 
the pole is preserved, and the theory is internally inconsistent because it is impossible to 
determine g(A) for all A (Ref. [IE]). In the latter case the position of the pole is given by 
the equation 




(24) 



dg 



(3(g) = (3 g 2 + Pig 3 + ■■■ , 



(25) 



din A 2 




(26) 
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which, for small values of g, leads to the result 



A 2 C = const • m 2 e 1//3 ° 9 , (27) 

which differs from (21) only by an insignificant constant factor. Thus, the existence of 
renormalon singularities seems fairly convincing for internally inconsistent theories. Con- 
versely, there is no reason for them in "good" theories. 

Since the behavior of the function (3(g) at g ^ 1 is unknown, the presence or absence 
of renormalon singularities is a matter of belief. We stress, however, the following point. 
Factorial contributions of individual diagrams exist in all field theories in which the ex- 
pansion of the j3 function (25) begins from the quadratic term: then the interaction on 
the scale is described by a formula of the type (20) with the replacement of A by k, 
whose expansion gives ((3olnk 2 ) N in the A^th order [see (13)]. To resolve the question 
of the internal inconsistency of a theory, one should know all the coefficients in the expan- 
sion (25). Therefore, it would be incorrect to consider the formal existence of renormalon 
contributions as an indication of the internal inconsistency of a theory. 



3. ANALYTIC PROPERTIES OF THE BOREL TRANSFORMS OF cp 4 

THEORY 



3.1. Expansion of the class of Borel transformations 

For the ensuing treatment it is convenient to expand the class of Borel transformations, 
setting 



F(g) = J dxe- x x h °- x B{gx), 



OO T? 



7 In particular, the result (27) is valid when the expansion (25) is truncated at a finite number of terms, 
provided the polynomial obtained is positive. On these grounds it is easy to draw the erroneous conclusion 
that the high-order terms of the expansion of the p function are insignificant. Parisi's arguments [161 117) 
regarding the momentum dependence of Borel transforms exactly follow this line of reasoning. In fact, the 
character of the solution of Parisi's equations [17] depends significantly on the behavior of (3(g) at infinity. 
In particular, they are easily solved for the model function (3{g) = flog 2 /(l + Xg) with A > 1 and lead to 
a result which differs qualitatively from its one-loop analog. 

8 The concrete sequences of the renormalon diagrams can differ somewhat in different theories. For 
example, in ip 4 theory the significant diagrams do not reduce to chains of "bubbles" (Fig. 2,c), but form a 
so-called parquet [35] , 
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with the arbitrary parameter 6 > 0, instead of (6) and (7). If B(g) and B(g) are the Borel 
transforms corresponding to the parameters b and b\ (for definiteness, we set b\ > b ), it 
is not difficult to derive the conversion formula: 

B(g) = ^ — . Jdx h \] B ( . (29) 

We define the analyticity region of B(g) by constructing the so-called Mittag-Leffier star 
[30J, i.e., by drawing cuts from all singular points to infinity along rays drawn through 
these points from the origin of coordinates. If g lies in the analyticity region of B(g), the 
integration path in (29) does not pass through its singularities and B(g) is also analytic. 
If g c is a singular point of B(g), the integration path in (29) for g = g c unavoidably 
passes through g c , generating a singularity in B(g). For the interesting case of power-law 
singularities we have the correspondence rules 

B(g) = AT{-0) {^jA - B{g) = AT(-/3 -b x + bo) ' " (30) 

for noninteger (3 + b\ — bo and 

B(g) = AT(-P) (i^lY^ B ig) = At£! ( hZl)\ (hZl) , (31) 
V 9c J n\ V 9c J V 9c J 

if j3 + b\ — bo = n is an integer. 

We see that the analyticity region for all the Borel transforms is identical and that it is 
sufficient to establish it for any fixed bo. The choice b = 1/2 is convenient for investigating 
functional integrals, since a simple result is obtained in that case for the Borel transform 
of an exponential function: 

F(g) = g-« - B{g) = °° S ^ = J_ {exp(2^ + c.c.}, (32) 

which preserves its exponential form. This permits writing an explicit expression for the 
Borel transform of the functional integral (1): 



Bi{g) = ^= J D(pexpi-S Q {(p}) 



exp 2iJgS int {(f} + c.c 



(33) 



The integrand is a regular function, and the analyticity region of Bj[g) is determined by 
the condition for convergence of the integral. 



3.2. Analyticity outside the negative semiaxis 
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For simplicity, let us consider scalar (p 4 theory. Generalization to the n-component case 
is trivial and reduces to only some complication of the notation. We assume that m 2 > 0, 
bearing in mind the subsequent analytic continuation to arbitrary complex m 2 . 

The integral (33) for ip A theory is defined well for positive values of g, since its conver- 
gence is determined by an exponential function of —So{(p} and is obvious after the Fourier 
transformation of <p(x): 



SoM = \jd d x {(V^) 2 + mV} = \ J2(k 2 



+ m )\ip k \ . 



(34) 



For the analytic continuation to complex g we turn the integration path in (33), setting 

g = get*, v = (pe-™l\ (35) 
where g and (p are real, and g > 0. Then the integral in (33) takes the form 

J D<p exp (-S {<p}e 



-i*/2 



exp 2iJgS int {(f} + c.c 



(36) 



and converges for — ir < ^ < n. Thus, the Borel transform is analytic outside the negative 
semiaxis. 



3.3. Analyticity within a circle 

We utilize the formal technique used in Refs. [3] and [H] and introduce the function 

SoW? 



i?M) 



Then (33) is rewritten in the form 

1 



Bi(g) 



20F 



D(p exp 



4S int M' 



1-i 



(37) 



1/2- 



c.c, 



and after the replacement of R{(f} by the constant Rq, it is analytic within the circle 
\g\ < Rq. Let us now suppose that 

RW} > Ro (39) 

for all (p, i.e., Rq is the exact lower bound of R{^}. Setting g = — \g\e n (— n < 7 < 71"), we 
have the inequality 



^)l<^/^(exp f- 



RW) 



1/2 



cos • 



7 



SoM +exp(-S M) , (40) 



which ensures convergence of the integral in (33) and, consequently, its analyticity within 
the circle \g\ < Rq. 
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To find -Ro) we consider the variational problem of minimizing R{tp}. It yields the 
equation 

-A<p(x) + m 2 <p(x) - dp\x) = 0. (41) 

where 

c = $)M 
2S- mt {ip}' 

which, after the replacement <p(x) — > ip(x)/y/C, transforms into the standard equation of 
an instanton of ip 4 theory. Using it, we can easily show that Ro = S{ip c }, which establishes 
the required analyticity region (Fig. 3,b). Questions concerning the absence of instantons 
in massive four-dimensional theory [UJ are discussed in Sec. 3.5. 

Apart from the integral (1), some other functional integrals containing products of the 
type <f(xi)<f(x2) ■ ■ ■ <p(xm) in the preexponential factor are of interest. The presence of 
such products does not influence the convergence, and all the proofs performed remain 
unchanged. 



3.4. Invariance relative to algebraic operations 

As 't Hooft pointed out [15] , the singularities of Borel transforms are not shifted when 
algebraic operations are performed on the original functions. This can easily be proved for 
a modified definition of the Borel transform (10), which differs from (6) and (28), since 

F(g) =F + F x g + F 2 g 2 + F 3 g 3 + . . . , B(z) = F 6(z) + ^ + + ^z 2 + . . . (42) 

and B(z) contains a 5-function singularity at zero. The transformation of (10) by means 
of the replacement g — > 1/z reduces to a Laplace transformation and allows inversion. It 
can be used to express the Borel transform of the product F 3 (g) = Fi(g)F 2 (g) in terms of 
the known Borel transforms of the factors: 

z 

B 3 {z) = J dz' B l {z')B 2 {z - z'). (43) 
o 

It can easily be seen that the 5-function singularity in B 3 (z) corresponds to the definition 
(42) and that the singular points for finite z coincide with the singular points of B\ (z) and 
B 2 (z) (see the analogous reasoning in Sec. 3.1). In particular, the Borel transform g n is 
the function z n-1 /r(n), which is analytic for integer values of n, and multiplication of the 
function by g n does not alter its analytic properties in the Borel plane. 
HF 2 (z) = l/F^z), then 

z 

5{z) = J dz' 'B 1 (z')B 2 (z - z') (44) 

Q 
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and the 5-function singularity on the left-hand side cancels out with the 5-function singu- 
larities in Bi(z) and B 2 (z). At finite values of z the right-hand side contains singularities 
corresponding to singular points of B\(z) and B 2 {z), which are absent on the left-hand side 
and, therefore, compensate one another. This is possible only if B 2 (z) has singularities at 
the same points as B\ (z). 

The proof of the analogous statements for linear operations, viz., summation, differen- 
tiation, integration, etc., is trivial. 

The standard definition of the Borel transform (6) is obtained from (10) and (42) when 
Fq = after the replacement F(g) — > gF(g). In this case the 5-function singularities 
disappear, and the remaining singularities are preserved at the same points due to the 
insignificance of the multiplier g. The definition (6) corresponds to the definition (28) with 
b = 1, and, by virtue of Sec. 3.1, the analysis performed can be extended to arbitrary b . 

Since all the quantities entering into the theory, viz. the Green's functions, vertex parts, 
etc., can be expressed in terms of functional integrals with identical analytic properties (see 
the end of Sec. 3.3) using a finite number of algebraic operations, their singular points in 
the Borel plane are the same as for the integral (1). 

3.5. Renormalization procedure 

The absence of ultraviolet divergences was implicitly assumed above. In ip A theory 
this is correct for d < 2. For 2 < d < 4 a continual theory without divergences can be 
constructed by introducing counterterms into the Lagrangian [4*ol I5T)] . In the simple case 
where only renormalization of the mass is required (2 < d < 4) the corresponding term in 
(4) is rewritten in the form 

ml(p 2 = (m 2 + Am 2 )y? 2 = (m 2 + Ag + Bg 2 + Cg 3 + . . .)y? 2 , (45) 

where the coefficients A, B,C, . . . are chosen so as to cancel the divergences. When coun- 
terterms are present, the analytic properties of integrals of the type (1) become more 
complicated, since the coupling constant appears not only in the combination g(p A , but also 
in the form of g(p 2 , g 2 tp 2 , etc. One of the directions of renormalon-related activity involved 
specifically the introduction of additional terms into the Lagrangian and tracing the renor- 
malon singularities appearing [TBI HH1 EH] . A question arises in regard to the cancellation of 
singularities in the case when the coefficients in front of the additional terms are selected 
so as to remove divergences, for which an unequivocal answer could not be obtained. 

A simpler route is to explicitly introduce regularization and to use renormalizat ion- 
group equations. Here we have in mind the so-called cutoff scheme ^T\: the vertex parts 
are calculated perturbatively as functions of the bare charge go an d the cutoff parameter 
A, then scaling functions which depend only on go are obtained, and, finally, renormalized 
vertices, which depend on the renormalized charge g, are constructed [50]. In this case the 
explicit introduction of counterterms is not required, but all the details associated with 
their presence are taken into account, since the fundamental possibility of eliminating the 
divergences is essentially used to write the renormalization-group equations. 
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The simplest way of regularization consists of substituting <pe(p)<p, where p is the mo- 
mentum operator, for the term (Vy?) 2 in (4), which is brought into the form —(pA(p = ^pp 2 ^p. 
If 

e(p) = e(-p), e(p)>0, (46) 

then both the entire structure of the instanton calculations [49J and the proofs presented 
above are preserved. The only change occurs in the equation of the instanton (41), which 
is brought into the form 

e(p)(p(x) + m 2 ip(x) - Lp 3 (x) = 0. (47) 

When the regularization 

e(p) =p 2 +p 6 /A 4 (48) 

is employed, the dependence of the action S{(f} on the instanton radius R in four- dimensional 
ip 4 theory has the form shown in Fig. 8. If A = oo, there is degeneracy with respect to 
the instanton radius in the massless theory [T], while there are no instantons in the massive 
theory [49J because of the monotonic dependence of S{ip} on R. At finite values of A a 
minimum appears on the plot of S{ip} versus R at m 2 > (the dashed curves in Fig. 8), 
and instantons appear in the massive theory. Their action S{ip c } determines the positions 
of the singularities in the Borel plane. For A — > oo and arbitrary m 2 > the value of S{(p c } 
tends to the instanton action of the massless theory, and the positions of the singularities 
do not depend on m.0 

The renormalization-group equations (in the Callan-Symanzik form) are valid for the 
vertices T L ' N with N free tails and L two-line insertions [50J: 

T L > N (g ,A) = 0. (49) 

Writing out three such equations with different L and N, we can express the scaling func- 
tions (3(go), vido), and ^(go) in terms of the vertices T L ' N (g , A) using algebraic operations, 
which do not shift the positions of the singularities in the Borel plane. In the limit A — > oo, 
where Eq. (49) is valid, the dependence of the scaling functions on A disappears [50], and 
their singularities in the Borel plane correspond to the massless theory. 

The Gell-Mann-Low function (3{go) defines the relationship between the renormalized 
charge g and the bare charge go . Let the functions F and F\ be such that F Q (g ) = Fi(g). 
The relationship between the corresponding Borel transforms Bq and B\ [in the sense of 
the definition (10)] can easily be found for an infinitesimal charge transformation, go = 

9 This dependence can easily be obtained by describing an instanton by two parameters, viz., its radius 
and amplitude, and performing a variational estimation of the action. In the theory of disordered systems 
this corresponds to the optimal-fluctuation method [35] . 

10 For dimensionalities 2 < d < 4 the influence of A on the properties of instantons is insignificant, and 
the role of renormalizations reduces to the fact that the instanton equation contains the renormalized mass 
[SB] , The dependence of S{ip c } on m is preserved in this case. 



d d ( N\ 
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S{<p c } 
A 




Aj < A 2 < A 3 < ... 



R 

Figure 8: Dependence of the action S{<p} on the instanton radius R in four-dimensional 
(/? 4 theory in the absence of regularization (solid curves) and for a finite cutoff parameter 
A (dashed lines). 
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g + 2(3(g)6A/A [see (25)]: 



B x {z) = B (z) + ^Jdy [B (y) + yB' (y))Bp{z - y), (50) 

o 

where Bp(z) is the Borel transform of the function j3(g)/g. Equation (50) is analogous 
to Eq. (43); therefore, the analytic properties do not change in the course of the charge 
transformation. 

The vertices T L,N diverge as A — > oo, but they become finite after separation of the 
divergent Z factors from them and the transition from the bare to the renormalized charge. 
Since the Z factors are, in turn, expressed in terms of the vertices T L,N (Ref. [SO]), the 
renormalized vertices have the required analytic properties. 

The dependence of the scaling functions on the renormalization scheme is given by the 
following conversion formulas [5T] : 

h^))=m dq{q) 



v(q(9)) = v(g) -13(g) 



dg 

dlnp(q) 



dg 



Ud(g)) = V2(g) - m dlnv / q) ■ (51) 

dg 

The conversion functions q(g), p{g), and p%{g) for standard renormalization schemes (sub- 
traction, cutoff etc.) are expressed in terms of the vertices T L)N . Therefore, the analytic 
properties of the scaling functions are identical in all the schemes. In the general case the 
analytic properties of the conversion functions require additional investigation. 



4. CONCLUSION 

The results in Sec. 3 rule out the existence of renormalon singularities in ip 4 theory. If 
the arguments in Sec. 2 regarding the relationship to a Landau pole are considered con- 
vincing, ip 4 theory cannot be internally inconsistent. The same conclusion can be drawn 
on the basis of solid-state applications: a reasonable model of a disordered system reduces 
exactly to ip 4 theory [361 139] , and the internal inconsistency of ip 4 theory would signify the 
impossibility, in principle, of obtaining a mathematical description of this model. There- 
fore, a revision of the results in Refs. [34] and [33], in which indications of the internal 
inconsistency of ip 4 theory were obtained on the basis of an approximate reconstruction of 
the Gell-Mann-Low function, is urgently needed. 

The results of Sec. 3 refer only to p 4 theory and cannot be extended directly to other 
field theories; however, along with the qualitative arguments in Sec. 2, they demonstrate 
the nonconstructiveness of the conception of renormalons as a whole. Therefore, it would 
be of interest to generalize the method of proof used in Sec. 3 to other cases. 
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In quantum chromodynamics (QCD) the renormalon doctrine presently prevails [2"U]- 
[30J. However, the specific details of QCD in this context have never been stressed. For 
example, 't Hooft [15], speaking about QCD, gives explanations within <p 4 theory and quan- 
tum electrodynamics. In recent publications [26j SO] the term "naive nonabelianization" 
appeared, which essentially means neglecting the specific details of QCD. On the other 
hand, in QCD there is a special reason for the belief in renormalons, which has a purely 
phenomenological character. The analysis of experimental data leads to the conclusion 
that the contribution of the high orders has a momentum dependence ocl/g 4 (Ref. [22]). 
This dependence can easily be obtained from renormalon graphs, but, as is generally as- 
sumed, it cannot be obtained within the instanton method. The latter is based on the 
results in Refs. [10] and [11] , according to which the instanton contribution is proportional 
to 1/g 18 . However, it can easily be seen that a term ocl/g 4 appeared in Refs. [TU] and 
|llj . but contained divergences which the authors found difficult to eliminate; this term 
was "transported" to the renormalon sector with the reasoning that it "contributes to the 
renormalon singularity . . .rather than to the instanton one" (Ref. [10J, p. 287). If there 
are no renormalon singularities, this contribution was simply discarded; therefore, no real 
calculation of the Lipatov asymptotics were made for QCD. 

This work was stimulated by lengthy discussions with P. G. Sil'vestrov, whom we thank 
for opposing the renormalon doctrine, his critical remarks, and general assistance in ac- 
quaintance with the situation. We also thank B. L. Ioffe, L. N. Lipatov, and the participants 
in the seminars at the Institute of Physical Problems, the P. N. Lebedev Physics Institute, 
the Institute of Theoretical and Experimental Physics, and the St. Petersburg Nuclear 
Physics Institute for their interest in this work and useful discussions. 

This work was carried out with financial support from the INTAS (Grant 96-0580) and 
the Russian Foundation for Basic Research (Project 96-02-19527). 



Appendix 

Construction of interpolation polynomial 

The polynomial of degree N which coincides with the function f(x) at the points 
x , xi, x 2 , ■ ■ ■ , xn, is defined by the Lagrange formula [53J: 

N ^ ip'(x k ) (x - x k ) ' 

ip(x) — (x — xo)(x — x\)(x — X2) ■ . ■ (x — xn), (-42) 
and the interpolation error is given by the expression 

R N (x) = f(x) - P N {x) = i—-^ ^{x), (A3) 
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Such divergences also appear in ip 4 theory, and a procedure for eliminating them is well known [3 



24 



where £ belongs to the interval (x ,xn). 

The function (18), which is of interest to us, behaves as q n at n< N with slowly varying 
q, so that In q ~ In N. Neglecting these slow variations and omitting the common multiplier 
in (16) and (18), we have 

/(*)=«• A„~^. (A4) 
Taking into account that \ip(x)\ < A N+1 in the interval < x < A, we obtain 

and the interpolation error is small for 

A < N/ In N. (A6) 

To investigate the dependence of the coefficient A N on the positions of the points Xk, we 
set ip(x) = Reip(x + iO), and calculating lnip(x + «0) using the Euler-MacLaurin formula, 
for ip(x) we obtain the expression 

*p(x) = F(x)sinG(x). (A7) 
In particular, for the power-law arrangement of points 

x k = (k/N) a A, k = 0,1,..., N, (A8) 

at « > 1 we have 

F(ar) = (-l) N y/x(A-x) exp{N [a{x/A) 1/a + In A - a] } , 

= nN(x/A)^ a . (A9) 
For a high-order coefficient of the polynomial (A.l) we obtain 

is 

(the sum is determined by the term with k — 1), and for a ~ In A" the coefficient 
can be factorially small or factorially large, depending on the relationship between a and 
In A, so that the required value of (A4) falls in the range of variation. Thus, the required 
polynomial (16) exists in the interval < n < n , where n ~ N/ In N . 
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